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O ■ ABSTRACT 

^ Aims. We study the geometry and kinematics of the circumstellar environment of the Be star k CMa in the Bry emission line and its nearby 
, . continuum. 

^ Methods. We use the VLTI/AMBER instrument operating in the K band which provides a spatial resolution of about 6 mas with a spectral 
resolution of 1500 to study the kinematics within the disk and to infer its rotation law. In order to obtain more kinematical constraints we also 
use an high spectral resolution Pa/S line profile obtain in December 2005 at the Observatorio do Pico do Dios, Brazil and we compile V/R line 
profile variations and spectral energy distribution data points from the literature. 
" " ' Results. Using diiferential visibilities and differential phases across the Bry line we detect an asymmetry in the disk. Moreover, we found 
that K CMa seems difficult to fit within the classical scenario for Be stars, illustrated recently by a Arae observations, i.e. a fast rotating B star 
close to its breakup velocity surrounded by a Keplerian circumstellar disk with an enhanced polar wind. Finally we discuss the possibility for 
K CMa to be a critical rotator with a Keplerian rotating disk and try to see if the detected asymmetry can be interpreted within the "one-armed" 
viscous disk framework. 



Key words. Techniques: high angular resolution - Techniques: interferometric - Stars: emission-line. Be - Stars: Keplerian rotation - Stars: 
individual (k CMa) - Stars: circumstellar matter 



1. Introduction 

The "Be phenomenon" is related to hot stars that have at least 
exhibited once Balmer lines in emission with infrared excess 
produced by free-free and free-bound processes in an extended 
circumstellar disk. There is now a strong evidence that the disk 
around the Be star a Arae is Keplerian (Meilland et al. l2006al) 
and that this dense equatorial disk is slowly expanding. On 
the other side there are also clear pieces of evidence for a 
polar enhanced wind. This was already predicted for almost 
critically rotating stars as for a large frac tion o f Be stars. 
Recently, Kervella & Domiciano de Souza (120061) showed an 
enhanced polar wind for the Be star Achernar whereas this 
Be star presents no hydrogen lines in strong emission. Thus, 
it seems that a significant polar wind may be present even if 
the star is still in a normal B phase, signifying this enhanced 
polar wind would not be related to the existence of a dense 
equatorial envelope. However many issues remain unsolved 
on the actual structure of the circumstellar envelopes in Be 
stars which probably depends on the dominant mass ejection 
mechanisms from the central star and on the way the ejected 
mass is redistributed in the nea r circumstellar environment. 
Recently Meilland et al. ('2006b') reported theoretical spectral 
energy distributions (SEDs), Bry line profiles and visibilities 
for two likely scenarii of the disk dissipation of active hot 
stars, and account for the transition from the Be to the B 
spectroscopic phase. 

K CMa (HD 50013, HR 2538) is one the brightest Be star 
of the southern hemisphere (V=3.8, K=3.6). It is classified 
as a B2IVe star, and the distance deduced from Hipparcos 
parallax is 230 ±30 pc. The measured vsini values range 
from 220 kms"> (Dachs et al. Il989t Mennickent et al. '2004; 
Okazaki 1997; Prinjai983) to 243 k ms'' (Zorec et a L2005i) , 
its radius is 6 Rq (Dachs et al. Il989t Prinia ll989l) and its mass 
is 10 M© (Prinjall989|). 
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We must mention that the mass and radius determination of a 
Be star is not an easy task. For instan ce if w e assume values 
of masses and radii from Harmanec d 1988 1) compilation, in 
agreement with Schaller et al. (1993) non-rotating evoluti onar y 
model s, for t he effective temperatures used by Popper (Il980l) . 
Prinja (Il989l) and Fremat 62005i) we obtain the table [T] 



inK 


Mass 
inMo 


Normal Radius 
inRo 


Radius from 
parallax in Rq 


20000 


6.60 


3.71 


7.25 (6.46 - 8.24) 


23100 


8.62 


4.28 


6.26 (5.59-7.13) 


25800 


10.72 


4.83 


5.59 (4.98 - 6.36) 



Table 1. Masses and radii determination for k CMa from 



Harmanec (Il988h compilation for th e effective t emper atures 
given by Popper ( fl98Q) . Prinja (Il989l) and Fremat ( l2005b . 



Thus, for a main sequence star the stellar radius should be 
smaller than the 6 R© we have adopted but on the other side, our 
radius estimate based on the parallax and the chosen V magni- 
tude from the correlation between th e brig htness and emission 
strength, as proposed by Harmanec (120001) . gives the range of 
radii comparable to our chosen 6 R© used in our modeling. 

The star exhibits a huge IR-excess and a strong emission in the 
hydrogen lines making a good candidate for the VLTI/ AMBER 
spectro-interferometer (Petrov et al. (12006 )) using medium 
spectral resolution (1500). Our aim is to study the geometry 
and kinematics of the circumstellar environment of this star as 
a function of wavelength, especially across the Bry emission 
line and to detect any signature of a possible asymmetry of its 
circumstellar disk as already observed through a violet to red 
pea ks rat io V/R ~ 1.3 by Dachs et al. (19921) and Slettebak et 
al. (Il992h . 
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2. Observations and data reduction 

Dedicated observations of k CMa were carried out during the 
night of the December, 26th 2004 with the three VLTI 8m 
ESO telescopes UT2, UT3 and UT4 (See Table[3]for the base- 
line configurations). The data were reduced using the amdlib 
(vl.l5)/ammyorick(v0.54) software package developed by the 
AMBER consortium. It uses a new data processing algorithm 
adapted to multiaxial recombination instruments called P2VM 
for Pixel To Visibility Matrix algorithm. The squared visibility 
estimator is computed from the basic observable coming from 
this algorithm that is the coherent flux (i.e. complex visibilities 
frame by frame multiplied by the flux) and the estimated fluxes 
from each telescope. The principles of the general AMBER 
data reduction are described in more details by Millour et al. 
(I2OO4.) and Tatulli et al. (2006). 

The complex coherent flux allows also to compute differen- 
tial phase, i.e. averaged instantaneous phase substracted from 
achromatic atmospheric OPD and a wavelength-averaged ref- 
erence phase. This means that the differential phase is the dif- 
ference between the phase of the source complex visibility and 
a mean OPD. This leads to an average differential phase equal 
to zero on the observed spectral window and the lost of the ob- 
ject's phase slope over the wavelengths. This technique allows 
however to retrieve partial information of the object's phase 
and is almost equal to the object's interferometric phase if we 
have some spectral channels where we know the object's phase 
is zero. 

Moreover, it also allows to compute "differential" visibility 
(as defined in Millour et al. 2006 and in Meilland et al. , 2006a) . 
i.e. the instantaneous modulus of the complex visibility divided 
by the averaged visibility on all the wavelengths excepted the 
work one. This leads to an average differential visibility equal 
to 1 in the continuum. It has the advantage over the "classical" 
visibility estimator to be almost insensitive to rapid frame to 
frame variations of visibility (due to vibrations or atmospheric 
jittering for example) and therefore one can expect the differ- 
ential visibility observable to be more precise than the classi- 
cal visibility estimator given the current vibrations in the VLTI 



infrastructure, and even though the continuum visibility infor- 
mation is lost in this observable. 

For more information, the differential data reduction is de- 
scribed in details in Millour et al. (l2005h and Millour et al. 
(l2006h . 

Reducing the k CMa data with a good accuracy was quite 
difficult to achieve. We encountered specific problems related 
to this data set. Therefore, in addition to the tools furnished 
by the default package, some specific processing was added to 
reach the best precision on the interferometric observables. 



Table 2. Calibration stars diameters estimate d from spectro- 
photometric indices (computed as in Bonneau et al. l2006h and their 
associated errors. 



Star 


Diameter (mas) 


Error (mas) 


HD 75063 


0.50 


0.08 


HD 93030 


0.454 


0.006 



- First of all, no specific data were available to calibrate the 
fringe contrast of k CMa. We therefore looked at calibration 
stars observed during the same night for other stars and cor- 
rected their visibilities averaged overall the [2.1 3-2.2 l]/im 
observed spectral range from their estimated diameters (see 
table |2|i in order to monitor the instrumental-!- atmospheric 
transfert function (see Fig.[Tll- This transfert function is the 
visibility of a point source measured by the instrument, al- 
lowing us to correct the raw visibilities on the science star 
from the instrumental-specific visibility loss. The scattering 
over the time of the visibilities gives the dispersion due to 
the instrumental drifts and atmospheric fluctuations during 
the observing time. This leads to a visibilities dispersion 
estimate of 0.05 for each star, which leads to an error on 
calibrated visibilities of 0.07 ( V2 x O.OS^). 
Then we interpolate the estimated transfert function to the 
time of the science star observations (as in Perrin et al. 



20031). By using this technique, we find that the [2.13-2.21] 



yum averaged visibility of k CMa is really close to 1.0 with 
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an uncertainty of 0.07 on all the observed base lengths. This 
would normally be killing for the wavelength-dependence 
study of the visibilities, but as explained before, we expect 
to have differential visibility and differential phase estima- 
tors much more precise than the visibility estimator 
The lack of dedicated calibration star for k CMa should 
have lead to an unfeasibility to spectrally calibrate the dif- 
ferential observables, but fortunately an other calibration 
star (HD93030) was observed two hours later within al- 
most the same spectral window, which means that the spec- 
trograph grating did not move but that the detector win- 
dow was not exactly the same as for k CMa, allowing us to 
use the intersecting spectral channels between the two ob- 
servations without any calibration problem. Detailed data 
analysis of calibration stars tends to demonstrate that the 
main pattern on differential observables comes from a fiber- 
injection pattern (i.e. a pure internal AMBER instrumental 
effect) and that it is stable over several minutes in the 10"^ 
range for the differential visibilities and 10"^ radians for the 
differential phase at medium s pectra l resolution (R ~ 1500, 
see for instance Vannier et al.. 



HD 75063 



K CMa 



HD 93030 



2005) 



This allowed us to calibrate correctly the differential visi- 
bility and the differential phases (See Fig.|2]). In order to ensure 
our calibration, we checked that all the features mentioned in 
this article are akeady present in the uncalibrated data, and not 
added by pure noise-effects produced by the calibrator star 

We could expect to see an effect in the closure phase, but 
its modulation seems to be of the order of amplitude of the 
error bars ( 3 °or 0.05 radians), which means that we do not see 
any detectable signal in the closure phase. This non-detection 
confirms the result on the visibility and the low amplitude of 
the modulation on the differential phases: the object is almost 
non-resolved or barely resolved by the interferometer on the 
considered baselines (80m maximum). 

What we see in the observed data is a decrease in the dif- 
ferential visibilities in two of the three baselines of the order of 
0.07, much larger than the error bars (0.02 for the differential 
visibilities). This can be explained by an envelope larger than 
the star, visible in the emission line. 
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Fig. 1. Raw absolute visibilities of calibration stars corrected from 
their angular diameters and averaged over the [2.13-2.21] yum window, 
allowing us to monitor the instrumental+atmospheric transfer! func- 
tion (points respectively around Ih in red and 7h in blue). We have 
overplotted for comparison the raw visibilities of k CMa (around 3h 
in green). The k CMa visibilities have obviously the same value as 
the instrumental+atmospheric transfer! function within the error bars, 
leading to a calibrated visibility of 1, i.e. a non resolved or very barely 
resolved object on all baselines. 

We observe also a modulation in the differential phase of 
the order of 5°(0.09 radians), also higher than the error bars 
( 2°or 0.03 radians). The modulation of the differential phase 
show a "sine arch" shape, typical of a rotating object or a bipo- 
lar outflow but also show an asymmetry, mainly on the baseline 
UT2-UT3 (Bl). 

In order to obtain more kinematical constraints the star has 
also been observed in the J2 band (1.2283-1.2937yum) with the 
1.6 m Perkin-Elmer telescope and Coude spectrograph (with 
R=10000) at the Observatdrio do Pico dos Dias, Laboratorio 
Nacional de Astroffsica (LNA), Itajuba, Brasil. The spectra 
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were recorded in the night 20/21 November 2005, at seven dif- 
ferent positions along the slit using the Camara Infravermelho 
(CamlV) detector The images of the dark have been subtracted 
from each star's spectral image, wavelength calibration image 
and five flat-field images. To get the sky image we made the 
median combination of the star's spectral images (divided pre- 
viously by the average of flat-field images). The sky image has 
been subtracted from stellar images and the one-dimensional 
spectra were extracted and calibrated in wavelength using the 
standard IRAfQ procedures. Finally the continuum normaliza- 
tion around the Pa^ line has been performed using our soft- 
ware. The average profile of the line, which has been used to 
constrain the kinematics within the disc, is plotted in Fig.|2]and 

13 



3. Envelope extension and flattening 

In this section we present the AMBER data in order to obtain 
an estimate of the k CMa's envelope geometry and extension. 
Assuming that the measured visibility in the continuum, Vc, is 
only due to the central star and its circumstellar envelope and 
that the envelope is optically thin in the continuum, we can 
write: 



(1) 



where V^.^ and F^.^. are respectively the envelope visibility and 
flux in the continuum, V* and F* are the star visibility and flux 
in the continuum and F^. - F^c+^-k- 

The total flux is normalized, i.e. F^. = F^.^, + F^, = 1. Since 
the star is fully unresolved (pi, ~ 0.25mas (assuming a 6 R© 
seen at 230 pc) which corresponds to V* > 0.99 for the longest 
baseline at 2.1 yum, we assume in the following that V* = l. In 
order to estimate V^c we still have to determine the star and 
the envelope contributions at 2.1 jim. Using the fit of the SED 
given in Fig. [3] we estimate that at this wavelength the stellar 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under cooperative agree- 
ment with the National Science Foundation. 



emission is similar to the envelope contribution, i.e. F* = F^^ - 
0.5. 

We have the same relation for the visibility in the Bry line: 



VerFer + VcF, 



(2) 



where and F^ are respectively the measured visibility and 
flux in the Bry line. and F^. are previously defined and V^r 
and Fp, are the visibility and flux only due to the Bry line, 
i.e. without any stellar contribution and envelope continuum. 
Using the AMBER Bry emission line profile plotted in Fig. |2] 
and neglecting the underlying broadened photospheric absorp- 
tion line, we obtain Ff,— 0.5 and F,.=1.5 at the center of the line. 

The corresponding visibilities, deduced from Eq. 1 and 2 
and from the measurements shown in Fig.|2] are given in Table 
[3] Using a uniform disk model for the envelope contribution, 
for each measurement, we also estimate in Table [3] the corre- 
sponding angular diameters in the continuum and in the Bry 
line. Since the envelope is marginally resolved in the contin- 
uum we simply put an upper limit to its angular size. 

The envelope extensions in Bry given in Table [3] are 
strongly dependent on the sky-plane baseline orientation as 
seen in Fig. |4] where we plotted the k CMa (unresolved star 
+ uniform disk) model diameters as a function of the baseline 
orientation. 

The K CMa circumstellar disk seems to be elongated along 
Bi but since we only have 3 visibility measurements we cannot 
accurately determine the angular position of the major-axis as- 
suming an elliptical circumstellar disk. The envelope flattening 
given by the semi-major and semi-minor axis ratio, is about 
2+0.7. Assuming that the disk is geometrically thin (i.e. its 
opening angle is only a few degree) we can estimate the range 
for the inclination angle i: 39°<i<68°. The lower limit of 39° 
relies on the lack of constraint on the disk opening angle. 

4. SIMECA modeling 

In order to obtain quantitative fundamental parameters of the 
central star and its circumstellar disk, we used the SIMECA 



code developed by Stee (1 1994 and Stee & Bittar (.200 li) to 
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Base n'' 


1 
1 


L 


-1 

5 


Baseline 


UT2-3 


UT3-4 


UT2-4 


Length (m) 


4Z. / 




on s 


RA. (°) 


jl.D 


Izo 


y 1 .V 


\7 

Vc 










0.93+0.02 


0.95±0.02 


0.9±0.02 




0.85<V,.<0.95 


0.87<V,.<0.97 


0.82<V,-<0.92 


V 


>0.86 


>0.86 


>0.86 


v„ 


0.69<V„<0.85 


0.75<V„.<0.91 


0.60<V,,<0.76 


(j)ec (mas) 


<3.6 


<2.6 


<1.9 


0„- (mas) 


3.7<0,r<5.5 


2.O<0„<3.6 


2.6<0„<3.4 


(pec (D*) 


<15.5 


<11.2 


<8.2 


(per (D*) 


15.9<0„<23.7 


8.6«^„<15.4 


11.2 <(^„-<14.6 



Table 3. Bry Visibilities measured in the continuum (V^) and visibil- 
ity drop within the Bry line (V^/V^.). V,- calculated from the measured 
Vc and V,/Vc ratio. Deduced envelope contribution in the continuum 
(y ec) and in the line (Vf,-) for each baseline. The corresponding angu- 
lar diameters in the Bry line ((p„) and the nearby continuum (0^^) are 
computed using a uniform disk model for each envelope measurement. 
The corresponding extension in stellar radii are also given, assuming 
a 6R0 star at 230pc. 



model the k CMa circumstellar environment. Since this code 
was axi-symmetric we made substantial modifications in order 
to introduce a longitudinal dependence of the envelope density 
as evidenced from the AMBER data plotted Fig. |2l To con- 
strain the kinematics within the disk we use a Pa/J line pro- 
file obtained in December 2005 at the Observatorio do Pico 
do Dios, Brazil and plotted in Fig. |5] This profile is strongly 
asymmetric with a V/R double peak ~1.3. This V/R > 1 is 
usually interpreted in terms of a viscous disk similar to accre- 
tion disks where the gas and angular momentum are diffused 
outward by magnetohydrodynamic viscosity (Lee et al. 1991). 
Considering th e time -dependent structure of isothermal viscous 
disk, Okazaki ( 119971) showed that "one-armed" density waves 
can propagate within the disk and should reproduce the ob- 
served V/R variations from V/R >1 to V/R<1 seen in the line 
profiles (Hummel & Hanuschik I1997I) . Such variations were 
detected for many Be stars, w ith pe riod from a few years to 
over a decade (Hanuschik et al l 19951: Telling et al. il994 ). But 



in the case of k CMa the V/R rat io has remained const ant for 
the last twenty years (Dachs et al. ll992t Slettebak [l992l) . 

In figure 5 we over-plotted the supposed "symmetric 
part" of the Pa/3 line profile, using an axi-symmetric model, 
and the asymmetric residual that may be produced within 
the "one-armed" oscillation over-density. This effect must be 
compatible with the asymmetric differential phase variation 
across the Bry line for the Bi baseline plotted in the bottom 
part of Fig.|5]since the emitting regions in Pa/? and Bry must be 
very close each together. The asymmetric contribution to the 
Bry emission is about 20 to 30% of the total emission in this 
line whereas the mean projected velocity of the inhomogeneity 
is -130±20 kms"' . Using a SIMECA model at 230pc we deter- 
mined that the projected separation between this over-density 
photocenter and the central star is about 6.5*. 

The parameters obtained for our best model are given in 
Table|4]with the coiTesponding spectroscopic and interferomet- 
ric observables of Fig. |2] This best model includes an over- 
density along the disk major axis at -h20°, coiTesponding to an 
over-luminosity of 30% of the total flux in the line, and the 
agreement with the VLTI/AMBER data, the SED (Fig. [3]l and 
the PayS line profile is very good as it can be seen Fig. |2] The 
very nice agreement with the differential visibility and phase 
across the Bry line for the three bases validates the chosen disk 
geometry and kinematics. The 2. ljum continuum visibilities ob- 
tained with the 3 baselines, respectively Vi=0.92, V2=0.96 and 
V3=0.94 are also compatible with the 0.93 lower limit mea- 
sured with AMBER. The corresponding continuum intensity 
map in the continuum at 2.15 /vm is plotted Fig.|6l The evalu- 
ation on the uncertainties of the parameters of our model is far 
below the scope of this letter and will be studied in depht as 
soon as we get more constraining data. 

5. Discussion 

Following recent VLTI/AMBER and VLTI/MIDI observations 
of a Arae Meilland et al. (l2006ah concluded that this classical 
Be star fits very well within the classical scenario for the "Be 
phenomenon", i.e. a fast rotating B star close to its breakup 
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parameter 


value 


Tts 


22 500K± 1000 


Radius 


6R0 ± 1 


Inclination angle i 


60° ± 10 


Equatorial rotation velocity 


240 km S-' ± 20 


rotation law exponent 


0.32 ± 0.1 


Photospheric density ipphot) 


4. 10""g cm"^ ± 2 10"" 


Equatorial terminal velocity 


1 km s"' ± 10 


Polar terminal velocity 


1000 km S-' ± 100 


Polar mass flux 


2 1O-"M0 year-' sr"' ± 0.5 10"" 


ml 


10 ±5 


ml 


10 ±2 


CI 


30 ± 10 


Envelope outer radius 


23 R* ± 2 


Major axis position 


+28° ± 5 


Over-density position 


along the disk major axis 



Table 4. Parameters for the k CMa central star and its circumstellar 
environment for the best axi-symetric model 



velocity surrounded by a Keplerian circumstellar disk with an 
enhanced polar wind. This scenario was also confirmed for th e 
Be star Achernar by Kervella & Domiciano de Souza 
using VLTI/VINCI data, even if, for this latter case, the star 
was not in its active Be phase, i.e without any strong emission 
Une and no circumstellar disk. Nevertheless, Achernar was still 
a nearly critical rotator and was still exhibiting an enhanced 
polar stellar wind. We will see in the following that k CMa 
does not fit very well within this classical scenario. 

5.1. ISK CMa a critical rotator ? 

If K CMa was rotating close to its breakup velocity, i.e. 
Vc=463kms"', the inclination angle would be around 28° 
in order to obtain a measured vsini =220 kms"'). With this 
inclination angle the maximum flattening corresponding to 
a geometrically very thin disk is 1.12. Since we measure a 
flattening of about 2+0.7 this inclination angle can be ruled 
out. In our best SIMECA model the star is rotating at only 
52% of its critical velocity. We may argue that the measured 



elongation is not the envelope major axis but rather the 
enhanced polar wind. In this case the projected axis of the 
Be envelope is not identical to the rotation axis of the star 
Nevertheless, in order to obtain an asymmetry in the jet we 
need an extended optically thick disk, perpendicular to the 
jets directions, that may screen a least one part of the jet-like 
structure. Such extended optically thick disk should have been 
detected with the AMBER instrument which is not the case in 
our data. 

The value of the projected rotational velocity for an early-B 
star can be systematically affected by pseudo-photosphere, 
unrecognized optica lly thick parts of the Be envelope as shown 
by Harmanec ( 20021) for y Cas. He obtains for this star a vsini 
of 380 kms ' ins tead o f the often quoted value of 230 kms ' 
from Sletteback (Il992h . Nevertheless, takin g the l argest value 
for K CMa from the literature from Zorec (l2005h who found 
a vsini=243kms"' we still obtain an inclination angle of 32° 
which again is not in agreement with our measured flattening. 
Of course, if the discrepancy between the measured vsini 
and the "real" one is larger it may be possible that k CMa 
is still a critical rotator but it requires a factor of 2 between 
the measured and the tr ue vsin i which we found unrealistic. 
Even if Townsend et al. (2004) include the gravity darkening 
effect on the vsini values of rigid early-type rotators, assuming 
a rotation rate Q/D.^ of 0.95, they conclude that classic vsini 
determinations for BO to B9-type stars can be underestimated 
by 12 to 33%, f ar froin a factor of 2. Moreover, a recent paper 
by Fremat et al. (120051) studying the eff'ect of the gravitational 
darkening on the determination of fundamental parameters 
in fast rotating B-type stars found that on average the rate of 
angular velocity of Be stars attains only Q/Qc~0.88. 



In this last paper, Fremat et al. (l2005h estimate k CMa's effec- 
tive temperature to be 25790 + 713 K, a value significantly 
larger than the 2 2500 K used in our modeling. Moreover, 
Harmanec (l2000h found a positive correlation between the 
emission strength and brightness in the optical. Therefore we 
may use the minimum observed V magnitude of about 3.5 
to estimate the radius of the central star Combining with the 
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Hipparcos parallax and its error we obtain a radius between 
9 and 14 solar radii. Using the T^y^y of 25790 K and a radius 
of 14 Rq we obtain a stellar luminosity larger by a factor 8 
compared to our modeling and thus it is not possible to obtain 
a good agreement with the observed SED plotted Fig|3] Thus 
we are more confident with our 6 Rq used for our modeling 
and our finding that k CMa seems not to be a critical rotator. 
Nevertheless, regarding the uncertainties and the large errors 
of all measurements the breakup velocity cannot be totally ex- 
cluded. 



5.2. Is the rotation law within the disk Keplerian ? 

A Keplerian rotating law would produce a narrower double- 
peak separation in the Pa/3 line profile. Using a simple 
axi-symmetric Keplerian disk model the double-peak sep- 
aration would be about 90kms ' whereas we measure an 
asymmetric double-peak separation of about 160kms"'. Even 
if we subtract the emission of the over-density producing a 
larger double-peak separation by contributing to the V peak of 
the emitting Pa/? line, we still obtain a double-peak separation 
of about 120kms"' (see Fig.|5]l. The exponent of the rotation 
law used for our best SIMECA model is 0.32 whereas it should 
be 0.5 for a purely Keplerian disk. 

We may argue that Be stars vary strongly in time and thus line 
profiles shapes are time dependent. For instance actual Ha line 
profiles show a str ong e mission with a single peak whereas 
Bahng & Hendry (1 19751) saw a double-peaked Ha emission 



line, with the same double-peak separation of 160 km"' we 
obtained for Pa/3, with a shell core on their high-dispersion 
spectra. Nevertheless, these line variations are certainly related 
to the formation and disappearance of the circumstellar disk 
aroun d the sta r as shown by Rivinius et al. (i2001 ) and Meilland 
et al. (l2006bl) . Whatever the model is, a double-peak line pro- 
file is a clear signature of an extended rotating disk, at least 
if the kinematics is not dominated by a strong stellar wind i n 
the equatorial region as shown by Stee & de Arailjo (11994 . 
This double-peak separation i s a go od indication of the disk 
extension as shown by Huang (Il972h : Hirata & Kogure (,1984 ) 



and Stee & de Arailjo ( 19941) . We measure Vdisk sini at the disk 
outer radius (Rdisk) from the peaks separation, where v^,,,, is the 
rotational disk velocity at Rdisk- Thus we can write: 

^ Rtlisk \ ^ 



Vdisk sin(0 = V* sin(0 



(3) 



where v* is the star rotation at its photosphere. 



Assuming a Keplerian rotation (/3 -0.5) we obtain, using Eq. 3, 
R£//.s/(:=13.5 R* which is about 2cr from the 19.8 R* interfero- 
metric measurement, assuming that the measured elongation is 
the envelope major axis and not an enhanced polar wind (see 
discussion in the previous point). Note that these 19.8 R* found 
are obtained assuming a uniform disk for the envelope and thus 
is certainly a lower limit to the "true" disk extension in the Pa^S 
line. Thus it seems difficult to maintain a Keplerian rotation 
within the disk of k CMa. 



5.3. Is the "one-armed" viscous d/sk model a possible 
scenario for k CMa ? 

The asymmetry presently detected in the disk of k CMa 
seems to be hardly explained within the "one-armed" vis- 
cous disk framework. Following the viscous disk models 
by Okazaki ('1997I) and the observational detection of "one- 
armed" oscillations in the disk of ^ Tau by Vakili et al. (Il998h 
and y Cas by Berio et al. (Il999l) . the precessing period (P) of 
such oscillations should be confined within a few years up to 
about twenty years for the longer ones. On the other side, we 
tried to compile all the observational data available in order 
to obtain a "quasi-period" for the V/R variations. First we 
must mention that the V/R variations occur during the time 
intervals of observable presence of Be envelopes and that they 
can sh ow lon g-term, medium-term as well as rapid changes 
(Dachs ll98lh . Moreover, the very strong Ha line profile is not 
really suitable to V/R measurements since it is single-peaked 
and the illusion of apparent V/R changes can be related to the 
presence of telluric lines. Finally, compiling the data b etwee n 
1965 and 2003 for k CMa respectively from Jaschek d 1965b . 
Slettebak ( 1982) . Banerjee (i2000i) and this work, we were not 
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able to deduce an estimation of a quasi-period (Fig. IT). Several 
authors suggested a very long V/R variation (i.e. Okazaki 
P>28 years). This is a possibility but an equally plausible 
possibility is that the star had two episodes of V/R changes 
with much shorter cycle length separated by a period of 
quiescence documented by Dachs et al. il992t Slettebak 1992. 
More observations are needed since, if this first possiblity 
could be confirmed, then there is a problem for the one-armed 
model. This "pseudo-period" would be too long compared to 
theoretical predictions which can be hardly longer than two 
decades for a disk with a radius ~ 23R* (Okazaki, private 
communication). Moreover, the fact that this over-density 
remains confined along the major axis of the disk seems to be 
a very nice coincidence... 

Clearly, further observations are mandatory to confirm these 
conclusions or to see if other physical phenomena occurred 
within the circumstellar disk of k CMa. 
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Fig. 3. K CMa's Spectral Energy Distribution (SED) from SIMBAD 
CDS (triangles). Dotted line: emission from the central star assuming a 
black body with i?* = 6i?o, T^//=22500K and d=230 pc. Dashed line 
: free-free and free-bound envelope contribution from the SIMECA 
code between 0.3 and 100 jum . Plain line : Central star emission -i- 
envelope contribution. 
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Fig. 2. From the top to bottom : Pa/8 line profile from the 
Observatorio do Pico dos Dias, Brazil (dotted line) with our best 
model fit (plain line), Bry line profile, diff^erential visibilities and dif- 
ferential phases for the three baselines. For each plot, the dots with 
errors bars are VLTI/ AMBER data and the solid line is from our best 
SIMECA model (see section 4). 
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Fig. 4. K CMa diameters in the Bry line, assuming an vmresolved star 

+ uniform disk models, as a function of the baseline position angle (in 
mas). The length of each plot corresponds to the error bar measure- 
ment whereas diameters are given by the center of each error bar. 





Fig. 5. Upper picture: k CMa Pa/S line profile observed in December 
2005 at the Observatorio do Picos dos Dias, Brazil (solid line). 
Estimated symmetric part of the Pa/? profile (dotted line) using an 
axi-symmetric model. The asymmetric residue corresponds to the 
emission of "one-armed" over-density (dashed line). Bottom picture: 
Differential phase variation measured along the Bi baseline (dots with 
errors bars) and theoretical phase from the SIMECA code. 
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Fig. 6. Intensity map in the continuum at 2.15 jim. obtained with 
SIMECA for our best model parameters. The inclination angle is 60°, 
the central black dot represents the k CMa photosphere (0.25 mas), the 
bright part in the equatorial disk is produced by the over-density which 
is oriented along the Bi baseline. This over-density is also responsible 
for a 30% emission excess in the asymmetric V part of the Bry line. 
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